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ABSTRACT 

We present and analyze a new M' detection of the young exoplanet (3 Pictoris b from 2008 
VLT/NaCo data at a separation of « 4 AU and a high signal-to-noise rereduction of L' data taken in 
December 2009. Based on our orbital analysis, the planet’s orbit is viewed almost perfectly edge-on 
( i ~ 89 degrees) and has a Saturn-like semimajor axis of 9.50 AU/ 3 '®^^/. Intriguingly, the planet’s 
orbit is aligned with the major axis of the outer disk (0 ~ 31 degrees) but probably misaligned with 
the warp/inclined disk at 80 AU often cited as a signpost for the planet’s existence. Our results mo- 
tivate new studies to clarify how 8 Pic b sculpts debris disk structures and whether a second planet 
is required to explain the warp/inclined disk. 

Subject headings: stars: early-type, planetary systems, stars: individual 8 Pictoris 


1. INTRODUCTION 

Two decades of studies have argued that the nearby, 
12 Myr-old A- type star 8 Pictoris likely harbors a young 
planetary system (e.g. Smith and Terrile 1984; Kalas and 
Jewitt 1995; Mouillet et al. 1997); recently, Lagrange et 
al. (2009a, 2010) detected a « 9 ± 3 Mj planet around 
this star (/? Pic b). Imaged at projected separations of 
~ 6 AU and ~ 8 AU (November 2003 and October- 
December 2009, respectively), 8 Pic b - along with HR 
8799e (Marois et al. 2010; Currie et al. 2011) - may also 
provide a more direct comparison to the solar system’s 
gas giants than other directly imaged planets which are 
at wider separations (e.g. Fomalhaut b and HR 8799bcd 
Kalas et al. 2008; Marois et al. 2008). Current studies 
have yet to detect the planet at projected separations 
< 0.3" (Lagrange et al. 2009b; Fitzgerald et al. 2009). 
Data at these smaller separations could provide crucial 
constraints on the planet’s orbit. 

Imaging planets at small, < 0.3" separations re- 
quires significantly reducing quasi-static speckle noise 
and wavefront errors induced by imperfect AO correc- 
tions. Advanced observing/image processing techniques 
like angular differential imaging (ADI) coupled with PSF 
subtraction from a locally optimized combination of im- 
ages (LOCI) algorithm (Marois et al. 2006; Lafreniere 
et al. 2007) significantly attenuate speckles and increase 
sensitivity. New instrumentation, such as the Gemini 
Planet Imager (GPI Macintosh et al. 2008), will achieve 
far superior wavefront control. Generally, these efforts 
focus on planet imaging in the near-IR. However, rnid- 
IR imaging naturally overcomes some of these challenges 
as the achievable Strehl ratio is better and the planet-to- 
st ar contrast is most favorable. High-contrast imaging 
with Strehl ratios >0.9 can yield at least some planet 
detections close to the telescope diffraction limit (e.g. 
2 A/D for HR 8799d; Serabyn et al. 2010). Because M 
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band imaging often achieves these high Strehl ratios (Mi- 
nowa et al. 2010; Hinz et al. 2006), it may be a promising 
route for detecting very young and self luminous plan- 
ets at small A/D separations, despite a much higher sky 
background. 

In this Letter, we report a detection of 8 Pic b at a 
separation of ~ 0.21" extracted from archival Mi band 
VLT/NaCo data taken in November 2008. We also 
present a high signal-to-noise L ' detection of 8 Pic b 
from December 2009 data first published in Lagrange et 
al. (2010). We combine these data with recent data from 
Bonnefoy et al. (2011) and Quanz et al. (2010) to better 
constrain the orbit and atmosphere of 8 Pic b. 

2. OBSERVATIONS AND DATA 

Our study originates from the need to test the 
ADI/LOCI reduction pipeline first presented in Currie 
et al. (2010) and updated in Currie et al. (2011) at sep- 
arations smaller than those where the pipeline had pre- 
viously extracted planet signals (r < 0.375"). Because 
8 Pic b’s reported projected separation in 2009 was « 
0.3" (Lagrange et al. 2010), we chose the now publicly 
available Lagrange et al. L' band data from December 
29, 2009 to test our code performance. Lagrange et al. 
(2010) discusses the details of the L' band observations. 
The total field rotation in units of the image FWHM was 
~ 3 A/D, sufficient for using our reduction pipeline. 

Figure 1 (top-left panel) shows our processed L' band 
image using the LOCI algorithm in annular regions of 
250 x FWHM (Na — 250) with reference images selected 
from frames with at least 0.5xFWHM field rotation (S = 
0.5). The planet is easily detected and is well separated 
from residual speckle noise. The planet signal-to-noise, 
determined from the dispersion in pixel intensity values 
in concentric annuli, is SNR ~ 21: about a factor of 4-5 
greater than from Lagrange et al. (2010). 

Motivated by this success, we searched for additional 
8 Pic data in the VLT /NaCo archive taken in ADI mode 
between 2003 and 2009, finding a set taken on Novem- 
ber 11, 2008 with the L27 camera. Most of these data 
were taken in sparse aperture masking mode in K s , Li 
and M' bands over a span of ~ 4 hours: these data were 
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mentioned in Lagrange et al. (2009b) as not providing 
•good constraints on the companion. However, we found 
~ 13 minutes of the M' data taken in ADI mode with- 
out aperture masking at various times in between the 
masking data. Over the course of the entire observing 
sequence, the parallactic angle changed by ~ 100 de- 
grees, or ~ 2.4-3 A/D at 0.2"-0.25": sufficient for image 
processing with our pipeline. 

Basic image processing of the M' band data followed 
steps outlined in Currie et al. (2011). After registering 
each image and subtracting off the smooth seeing halo, 
we Fourier filtered the data to remove residual low spa- 
tial frequency noise and masked any hitherto unidentified 
bad pixels previously lost in the seeing halo. We explored 
a range of LOCI parameter space, varying 6, Aft, and the 
ratio of the radial to azimuthal lengths of the subtrac- 
tion annulus ( g ). Because 3 Pic b is very luminous in 
the mid-IR (e.g. A L' ss 7.7; Lagrange et al. 2010, and 
Section 3 of this work) , we focused on “aggressive” LOCI 
settings of 8 = 0.25-0.5, Na = 200-300, and g = 0.3-1, 
which better remove residual speckle noise. 

Figure 1 (top-right panel) shows our best reduced M' 
image. 3 Pic b is clearly detected in the southwest quad- 
rant ~ 0.2-0.25" from the star (SNR ~ 6). Manually 
inspecting each image between the radial profile subtrac- 
tion and final image combination steps and examining a 
signal-to-noise map of the median-combined image also 
shows that the peak does not result from latent image 
artifacts. Slightly different settings for 5, Na, and g also 
yield significant detections (bottom panels). 

3. ANALYSIS 

Our new M' band detection and high signal-to-noise 
L' band detection allow new constraints on the planet’s 
orbit. To derive precise astrometry needed to investi- 
gate the planet’s orbit, we adopt the NaCo plate scale 
and orientation for the L27 camera from Bergfors et al. 
(2011): 27.1 mas/pixel and a north position angle of -0.6 
degrees. These values are nearly identical to those for the 
L27 camera quoted by Lagrange et al. (2009a) for 2003 
NaCo data and for the S27 camera from Ehrenreich et al. 
(2010) calibrated from Trapezium data acquired closest 
in time to the 3 Pic data: our astrometric results does 
not leverage on which calibration we use. 

To fine tune our measurements, we correct for the pho- 
tometric and astrometric bias induced by LOCI process- 
ing by comparing the imputed fluxes and positions of 
fake point sources added to registered images with com- 
puted fluxes and centroid positions obtained after LOCI 
processing (e.g. Lafreniere et al. 2007; Thalmann et al. 
2009; Currie et al. 2011). While we lack unsaturated 
data from this run to directly confirm the PSF shape, 
unsaturated M' data taken in prior runs such as that- 
for HD 158882 (March 2007) show that the AO-assisted 
NaCo M' PSF core is axisymmetric and well reproduced 
by a simple Gaussian intensity distribution. For the V 
band data, the astrometric bias is minimal, whereas 0 
Pic b’s measured radial separation in M' band is biased 
by about +0.5 pixels (0.013”). The position angle offsets 
for both data were minimal. 

We determine the \V band position to be at a separa- 
tion of r — 0.210 ± 0.027" and position angle of 211.49 
±1.9 degrees. The L' band position is at 0.326 ± 0.013" 
and 210.64 ± 1.2 degrees (Table 1). Here we conserva- 


tively assume an uncertainty in radial separation of one 
pixel for M' band and 0.5 pixels for the (higher signal- 
to-noise) V band data. The position angle uncertainty - 
determined from the dispersion in values using different 
centroiding estimates (e.g. cntrd.pro vs. gcntrd.pro) is 
0.25 pixels (0.7 masxr), or 1.2 and 1.9 degrees for L' and 
M', comparable to uncertainties for 3 Pic b by Lagrange 
et al. (2009a, 2010) and Bonnefoy et al. (2011). Assum- 
ing that 0 Pic is 19.3 pc distant (Crifo et al. 1997), the 
planet was at a projected separation of 4.05 ± 0.50 AU 
on November 11, 2008 and 6.29 ± 0.25 AU on December 
29, 2009. 

To determine the range of allowable orbits for 3 Pic b, 
we follow the method described in Janson et al. (2011) 
used to model the orbit of the low-mass brown dwarf 
companion GJ 758 B (Thalmann et al. 2009; Currie et 
al. 2010), somewhat similar to earlier analyses for 8 Pic 
b (Lagrange et al. 2009b; Fitzgerald et al. 2009). In this 
approach, we perform a Monte Carlo simulation com- 
paring the astrometry to predictions from randomly se- 
lected orbits, where we allow all orbital parameters to 
vary. The minimum x 2 value in our simulation is x 2 ~ 
1.23. Given our data’s weak constraints on the orbital 
acceleration and the degeneracies due to the unknown 
line-of-sight components of planet position and velocity, 
no single ’best’ orbital solution emerges. Rather, the 
best-fitting solutions describe an extended, well- defined 
family of solutions that all match the data equally well. 
We choose a cut of x 2 < 2-23 (x 2 < x 2 min + 1) to repre- 
sent the family of best-fitting orbits. We also consider the 
results for a cut of x 2 < 8 (x 2 for a 1-a deviation in each 
of the data’s degrees of freedom): a family of ’average- 
fitting’ solutions. From the set of models satisifying this 
criterion, we determine the median value of each param- 
eter, weighted by the ratio of the mean to current or- 
bital velocity for the corresponding orbit, and identify 
the weighted 68% confidence interval about the median. 
We include astrometry from the highest signal-to-noise 
data separated in time by more than ~ 3 months (Table 
!)• 

Figures 2 and 3 displays our Monte Carlo simulation 
results. For a y 2 < 2.23 cutoff (Figure 2), the range of 
best-fit orbital parameters (weighted median, [weighted 
68% confidence interval]) include a p = 10.99 [8.18, 15.88] 
AU, i = 89.47 [89.19, 89.69] degrees, e = 0.12 [0.03,0.31], 
and f l = 30.89 [30.57, 31.17] degrees. For nearly circular 
orbits (e < 0.1), the range in semimajor axes is much 
narrower ( a p ~ 8-12 AU). If we relax our fitting criteria 
to accept models with x 2 < 8 (Figure 3), we find a p = 
9.50 [7.80. 13.43 ] AU and e = 0.10 [0.02,0.23]. More 
importantly, the inclination and longitude of ascending 
node are still nearly single valued: i = 88.93 [88.06, 89.40] 
and fl = 31.32 [30.56, 32.12]. Thus, 3 Pic b’s orbit must 
be viewed almost perfectly edge on, consistent with that 
for 3 Pic’s debris disks, with a north position angle of ~ 
30.8 degrees for the outer debris disk (Kalas and Jewitt 
1995, see also Boccaletti et al. 2009) but inconsistent 
with the inner disk position angle, which is offset by ~ 5 
degrees (Heap et al. 2000; Golimowski et al. 2006). 

Unfortunately, the M' band observations of 3 Pic were 
taken with the star saturated within ~ 3 pixels (~ 0.6 
FWHM) of the centroid and there were likewise no un- 
saturated standard star observations. We derive a very 
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crude magnitude estimate by scaling the M' PSF of HD 
158882 to the unsaturated portion of the 8 Pic PSF and 
use HD 158882’s known brightness (K a = 5.09; K s -M' 
~ 0) to calibrate 8 Pic b’s brightness. We estimate A M' 
« 8.02 ± 0.50 (Mm- ~ 9.99), where we consider the un- 
certainties in our PSF fitting scaling, the dispersion in 
individual planet magnitude estimates drawn from sep- 
arate reductions, and the intrinsic signal-to-noise of our 
detection. We determine an IJ contrast of A L' = 7.71 ± 
0.06. Combining the V measurement with the K s band 
and [4.05] data from Bonnefoy et al. (2011) and Quanz et 
al. (2010), we have three good quality photometric points 
to investigate the family of possible solutions for 8 Pic 
b’s atmospheric properties. 

Figure 4 compares the 8 Pic b photometry to best- 
fit spectra for models with log (g) = 3. 5/4-4. 5 and T e ff 
= 1000-1800 K for a range of cloud prescriptions: the 
Model A and AE thick cloud prescriptions respectively 
from Currie et al. (2011) and Madhusudhan et al. (2011) 
that best fit the HR 8799 planet SEDs, the Model E cloud 
deck prescription appropriate for brown dwarfs (Burrows . 
et al. 2006), and a cloudless atmosphere. The values 
for these models are, respectively, x 2 = 24.8, 12.3, 20.4, 
and 43.2 for Models A, AE, E, and the cloudless case. 
The AE thick cloud model provides the best fit. Thick 
cloud models also produce redder L'-M' colors at high 
temperatures, similar to that estimated here (« -0.22 ± 
0.50), though our lack of a reliable M' photometric cali- 
bration precludes strong conclusions. Good photometry 
is available in only three filters, so we cannot yet say that 
8 Pic b has thick clouds like the HR 8799 planets (Currie 
et al. 2011; Madhusudhan et al. 2011). 

The range of gravities and effective temperatures are 
log(g) = 3. 5-4. 5 and T e ff = 1400-1800 K. The implied 
masses for these models range between 4.1 Mj and 19.2 
Mj and ages range between 1 and 27 Myr, broadly con- 
sistent with the planet mass (9 ± 3 Mj; Lagrange et 
al. 2010), stellar age (12 Myr; Zuckerman et al. 2001), 
and likely formation timescale (< 3-5 Myr; Currie et al. 
2009). Planet fluxes in the near-IR (1-1.65 ^m) and the 
3-3.5 fim methane absorption trough are highly sensitive 
to cloud structure (e.g. Currie et al. 2011; Madhusudhan 
et al. 2011). Thus, J or H broadband data and/or nar- 
rowband 3-3.5 jim data will be critical in breaking model 
fitting degeneracies. 

4. DISCUSSION 

We present a new detection of 8 Pic b in M' band at a 
separation of ~ 0.21" (a pro j ec ted = 4.05 AU) from archival 
VLT/NaCo data taken in November 2008 and a high 
SNR rereduction of L' data first reported by Lagrange 
et al. (2010), using these data to constrain the planet’s 
orbit and atmospheric properties. For orbits whose fit to 
the data yield y 2 < 2.23, we find that the 8 Pic planet 
has a semimajor axis of a p = 10.99 AUl^’iitu an d a 
moderate/low eccentricity (e < 0.31). Admitting orbital 
solutions with y 2 < 8, the parameter ranges are a p — 
9.50 AU/f^ ; u and e < 0.23. In both cases, values for 
the planet’s inclination (i ~ 88.06-89.69 degrees) and 
longitude of ascending node (fi ~ 30.56-32.12 degrees) 
are tightly constrained and imply that the planet’s orbit 
is almost perfectly aligned with the outer debris disk, 
but not the inclined inner disk (Q ~ 35-36 degrees). We 


cannot extract reliable photometry from our M' band 
data; new data at 1-1.65 /cm and ~ 3-3.5 gm is needed 
to constrain 8 Pic b’s atmosphere. 

Numerous studies of the 8 Pic debris disk(s) have iden- 
tified the star as harboring a young planetary system 
(e.g. Smith and Terrile 1984; Kalas and Jewitt 1995; 
Mouillet et al. 1997; Weinberger et al. 2003). More re- 
cently, the presence of a warp in the disk at ~ 80 AU - 
due to the combined effects of the main disk with PA ~ 
30.8 and a second, inclined disk offset by 5 degrees - was 
identified as a clear signpost of a perturbing planet (e.g. 
Mouillet et al. 1997; Heap et al. 2000; Golimowski et al. 
2006), motivating high-contrast imaging studies to im- 
age the planet. Lagrange et al. (2009a) identified 8 Pic 
b as a likely source of the inclined disk and used the disk 
morphology to derive mass estimates (see also Lagrange 
et al. 2010). 

Our results suggest that 8 Pic b is probably not aligned 
with the inner disk/warp but rather the main disk, as 
the allowed range in f l is offset from the main disk as 
measured by Kalas and Jewitt (1995) by no more than 
~ 1 degree. Furthermore, the planet may be misaligned 
with the submm disk emission (Wilner et al. 2011), which 
is sensitive to dynamical sculpting by planets (Kuchner 
and Stark 2010). However, models accounting for the in- 
clined inner disk presume that the planet’s orbit is also 
inclined relative to the main disk (e.g. Mouillet et al. 
1997; Augereau et al. 2001). New beta Pic b astrom- 
etry, a more precise astrometric calibration of existing 
beta Pic NaCo data by determining and correcting for 
image distortion, and a detailed relative calibration be- 
tween NaCo data and data revealing the disk will further 
clarify how beta Pic b’s orbital plane compares to that 
for the main disk and the inner disk/warp. Furthermore, 
our new orbital constraints for 8 Pic b strongly motivate 
new studies of the dynamical sculpting of 8 Pic’s debris 
disk by planets. If 8 Pic b or non-planet related mech- 
anisms (e.g. Armitage and Pringle 1997) fail to explain 
the inclined debris disk/ warp, the existence of additional 
planets in the system may be required. 

Our M' band detection demonstrates that it is possible 
to directly image planets at separations approaching the 
telescope diffraction limit without sparse aperture mask- 
ing interferometry (SAM; Ireland and Kraus 2008). The 
high Strehl ratio, large amount of field rotation, large 
mid-IR planet brightness, and LOCI processing pipeline 
are the keys to closing this gap. While SAM can de- 
tect planets interior to the telescope diffraction limit, it 
is overall less sensitive. However, the techniques can be 
complementary, yielding detections or robust limits on 
infant gas giant planets around the youngest stars on ~ 
5-100 AU scales. 

Upcoming facilities like GPI, SPHERE, SCExAO, and 
Project 1640 achieve higher contrast at small inner work- 
ing angles in the near-IR. primarily through more sophis- 
ticated wavefront control (Macintosh et al. 2008; Beuzit 
et al. 2008; Martinaehe and Guyon 2009; Hinkley et al. 
2011). Our results, coupled with previous L' band de- 
tections of 8 Pic from Lagrange et al. (2009a, 2010) and 
the high signal-to-noise L' band detection of HR 8799e 
(Marois et al. 2010) suggest that the mid-IR may also 
be fertile ground for new exoplanet detections at small 
separations for very young systems. Young, nearby 1.5- 
2 Mq stars like 8 Pic are particularly promising targets 



TABLE 1 

/? Pic Data Used in this Paper 


Astrometry 


Date 

Filter 

Separation ("), Position Angle (°) 

Reference 

11-10-2003 

V 

0.411 ± 0.008, 31.7 ± 1.3 

Lagrange et al. (2009a) 

11-11-2008 

M’ 

0.210 ± 0.027, 211.49 ± 1.9 

this work 

12-29-2009 

V 

0.326 ± 0.013, 210.64 ± 1.2 

this work 

03-20-2010 

K, 

0.345 ± 0.012, 209.8 ± 0.8 

Bonnefoy et al. (2011) 

Photometry 

Date 

Filter 

Absolute Magnitude 

Reference 

03-20-2010 

K s 

11.20 ± 0.12 

Bonnefoy et al. (2011) 

12-29-2009 

V 

9.73 ± 0.06 

this work 

04-03-2010 

[4.05] 

9.77 ± 0.23 

Quanz et al. (2010) 

11-11-2008 

« 9.99 ± 0.50 

this work 


Note. — Our astrometry and photometry are drawn from the three 
separate reductions shown in Figure 1. The L' measurement assumes 
L' g pic = 3.45. The M 1 photometry lacks reliable photometric calibration 
and thus is not useful for atmospheric modeling. 

for direct imaging surveys (e.g. Crepp and Johnson 2011) 
and many have resolved debris disks (e.g. HD 181327, 

Schneider et al. 2006; Chen et al. 2008). Imaging mas- 
sive planets in such systems can yield additional studies 
of planet-disk interactions, such as those motivated by 
this work. 
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Fig. 1. ~ (Top-left) Reduced VLT/NaCo L' band image from Dec. 29, 2009 data showing a 21-cr detection of j3 Pic b at a separation 
of ~ 0.32”. (Top-right) Reduced image showing the detection of fi Pic b in M' band (same image size). The LOCI parameters used to 
construct the reduced image include 6 = 0.25, g = 0.4, and Na = 250xFWHM. (Bottom panels) Reduced images using less aggressive 
LOCI settings - S = 0.5, <?=0.5, and Na — 300xFWHM (left panel) and S = 0.5, <7=1, and Na = 300xFWHM (right panel) - yielding 
detections of SNR ~ 5 and 4.5, respectively. In general, we detect the Q Pic planet at a 4-6 a level using a range of LOCI parameters: <5 
= 0.25-0.5, g = 0.3-1, N A = 200-300. 
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Fig. 2. Orbital analysis results following the Monte Carlo method outlined in Thalmann et al. (2009), showing parameters for solutions 
fulfilling x 2 ^ 2.23. The left panel displays the family of solutions in semimajor axis, eccentricity, and inclination space. The right panels 
show histogram distributions of these three parameters and the longitude of the ascending node. The vertical dashed line identifies the 
weighted median value for each parameter; the vertical dotted lines define the 68% confidence interval for each parameter. 
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Fig. 3. Same as previous figure except for orbital solutions fulfilling x 2 < 8. 



At/im) A (^m) 

Fig. 4. Figures comparing the 0 Pic b photometry to spectra for a range of cloud prescriptions: the Model A ’thick cloud’ limit used 
in Currie et al. (2011) (upper-left), the Model AE ’thick cloud’ prescription from Madhusudhan et al. (2011) (upper-right), the Model 
E brown dwarf cloud prescription from Burrow’s et al. (2006) (lower-left), and a cloud-free atmosphere also from Burrows et al. (2006) 
(lower-right). The blue dots show the flux predictions for each bandpass. All models assume a modal particle size of 60 y,rn. 










